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We propose a systematical approach to construct generic fractional quantum anomalous Hall 
(FQAH) states, which are generalizations of the fractional quantum Hall states to lattice models with 
zero net magnetic field and full lattice translation symmetry. Local and translationally invariant 
Hamiltonians can also be constructed, for which the proposed states are unique ground states. 
Our result demonstrates that generic chiral topologically ordered states can be realized in lattice 
models, without requiring magnetic translation symmetry and Landau level structure. We further 
generalize our approach to the time-reversal invariant analog of fractional quantum Hall states- 
fractional topological insulators, and provide the first explicit wavefunction description of fractional 
topological insulators in the absence of spin conservation. 
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Introduction- Topological states of matter arc quan- 
tum states which are distinguished from ordinary states 
by topological properties, rather than more conventional 
properties such as symmetry. The first examples of topo- 
logical states of matter discovered in nature are the inte- 
ger and fractional quantum Hall states P, Q which are 
realized in two-dimensional electron gas systems in a 
strong perpendicular magnetic field. The integer quan- 
tum Hall (IQH) state is characterized by a Chern number 
of the geometrical gauge field defined in the magnetic 
Brillouin zone[3[ or parameter space of twist boundary 
conditions 0, Hj]^ which determines the integer n in the 
quantized Hall conductance o~ xy — ne 2 /h. In a work in 
19880, F. D. M. Haldane proposed the first realization 
of the IQH state in a band insulator model without net 
magnetic field. During recent years, several semiconduc- 
tor systems have been proposed[7H9j which may realize 
such an IQH state without Landau levels, named as the 
quantum anomalous Hall (QAH) state. 

Based on the understanding of QAH states, a natu- 
ral question is whether one can also find fractional quan- 
tum anomalous Hall (FQAH) states, which are fractional 
quantum Hall (FQH) states without Landau levels. Dif- 
ferent from IQH states, FQH states necessary require 
electron-electron interaction, which makes the general- 
ization to the systems without Landau levels more dif- 
ficult. In ordinary FQH states the kinetic energy of 
electrons are quenched due to the flat Landau level, so 
that the electron interaction effect can be significant and 
lead to topological nontrivial states. On the contrary, 
in a QAH system the energy dispersion is in general not 
flat, so that if we consider the interaction effect, there is 
competition between kinetic energy and interaction en- 
ergy which usually disfavors the topological nontrivial 
states. Recently, specific QAH models with almost flat 
band dispersion and nontrivial Chern number has been 
constructed (lfj|. Numerical evidences of FQAH states 
have been found in such flat band models when interac- 
tion is considered [ill [l2| . However, the understanding of 
FQH states based on wavefunctions such as the Laughlin 



wavefunction [13| cannot directly apply to FQAH states, 
since the the single-particle and many-body wavefunc- 
tions in the QAH system are defined on lattice and can- 
not be written as analytic functions. 

In this Letter, we propose a systematical way to de- 
scribe the FQAH states by constructing model wavefunc- 
tions. We show that one-dimensional maximally localized 
Wannier functions can be defined in QAH states, which 
plays the same role as that of the Landau level wave- 
functions in IQH and FQH states. Based on the Wan- 
nier function basis, we construct the analogue of Laughlin 
wavefunctions in FQAH, and also obtain the analogue of 
the pseudo-potential Hamiltonians [13, [3 of which these 
wavefunctions are exact ground states. Once such a one- 
to-one mapping between Wannier functions in QAH and 
Landau level wavefunctions in QH is defined, each wave- 
function constructed in FQH has a counterpart in the 
FQAH. This demonstrates that the physics of FQH does 
not rely on any special property of the Landau level prob- 
lem, such as the wavefunctions being analytic functions, 
and the magnetic translation symmetry. Instead, frac- 
tionalizcd topological states exist generically in a flat (or 
nearly flat) band with a nontrivial Chern number. 

QAH states and one- dimensional Wannier functions. - 
The QAH state is a band insulator described by a Hamil- 
tonian 



H = J2<lh(k)c k 



(1) 



The single particle Hamiltonian h(k) is a N x N Her- 
mitian matrix for a system with N bands. Denote 
the eigenstatcs of the single particle Hamiltonian h(k) 
as \n, k) , n = 1,2, ...,N with eigenvalue E n , the Hall 
conductance of the system is determined by the first 
Chern number^ C\ = J d 2 \s.f xy (k) with f xy (k) = 



i) 



d y a x and a^(k) 



EE 



E n <0 



(n,k| di \n,k) , i 



x, y. The QAH state is described by a nontrivial 
Chern number. For a system with C\ = 0, the 
wavefunction of states |n, k) can be taken as single- 
valued in the Brillouin zone, from which one can con- 



2 



struct the Wannier function basis by a Fourier transform 
\W m ) = Ek,™ eik,x ^ m (fc)|m,k) with U nm (k) 



(a) 10r 



some smooth unitary transformation. On the contrary, 
for QAH states with C\ 7^ 0, it is well-known that Wan- 
nier function localized along both x and y directions can- 
not be defined since the wavefunction can not be taken as 
single- valued through the Brillouin zone.fli, 17 1 Instead, 
one-dimensional (ID) Wannier functions can be defined 
which are eigenstates of k y and are maximally localized 
along x direction [I3 - l20j . For each fixed k y , all states with 
momentum k y form a one-dimensional subsystem with 



the Hamiltonian Hn)(k y ) = Ea 



h[k x , fc y )cfe x fc 



For one-dimensional systems there is no obstruction in 
getting localized Wannier function. Maximally localized 
Wannier functions can be obtained as eigenstates of the 
projected position operator X = P_a;P_ with x the x- 
direction coordinate operator and P- the projection to 
occupied bands. 21 

For simplicity, we consider a QAH system with only 
one occupied band denoted by \k x ,k y ). The Berry's 
phase gauge field Oj = — i {k x , k y \ di \k x , k y ) with i = x,y. 
One can always make a gauge choice a y = 0, in which 
case the explicit form of the maximally localized Wan- 
nier function is 



\W(ky,x)) =i- 1 /2^ e -i^ ^(P*,ky)dp. 



i k j. 



(2) 



with 9{k y ) = Jq™ a x (px, k y )dp x - x € Z labeles the lat- 
tice sites. The phase factor with e lS ( k v) k ^/ 27T guarantees 
that the Bloch function is periodic for k x — >• k x + 2tt. 
In the gauge transformation \k x , k y ) — > e lL P( k *< k y) |/ Ca . j k y ), 
the Wannier function is gauge invariant up to an over- 
all phase: \W(k y ,x)) -» e 1 ^ ^ \W(k y ,x)). It can be 
verified directly that the center of mass position of the 
Wannier function \ W(k y ,x)) is given by 

(x) = (W{k v ,x)\x\W(k v ,x)) =x-e{k y )/2n 7 

so that 9(k y )/2-K is the shift of the Wannier function away 
from the lattice site, i. e. the charge polarization. [22l| 

Since the polarization 9(k y ) is determined by the Wil- 
son loop of the Berry's phase gauge field, the Chern num- 
ber in the Brillouin zone corresponds to the winding num- 
ber of 6[ky), i.e., C\ = — J Q % d6(k y ). For nontrivial 
Ci, the Wannier function \W(k y ,x)} is not periodic in 
ky, since its center-of-mass position shifts by C\ when k y 
is tuned continuously from to 27r. In other words, the 
Wannier functions satisfy the following twisted boundary 
condition: 

\W(k y + 2tt, x)) = \W(k y , x + C x )) (3) 
As an example, we consider the lattice Dirac model0,[23| 
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FIG. 1: (Colour online) (a) The center-of-mass position (x) of 
the Wannier functions versus k y . (b) (x) versus the extended 
wavevector K y defined in Eq. ([5]). (c)The profile of a Wannier 
function, (d) The profile of the coherent state wavefunction 
defined in Eq. (I10|) . All the results are calculated for the 
two-band model (TjlJ) with the parameters M = — 1,B — 1/2 
on a 50 x 50 lattice. 



which is a two-band model with the Hamiltonian 

h(k)= £ d a (k)a a (4) 

a=l,2,3 

with {d\, di, <^3)(k) = (sin k x , sin k yi M + 2B(2 — cosk x — 
cosfc y )). For B > and -2 < M/2B < 0, the system 
has Chern number C± = 1. The center-of-mass position 
of the Wannier functions can be obtained numerically as 
shown in Fig. Q] (a), which shows the shift of Wannier 
function position under k y — > k y + 2-k. As a consequence 
of this shift, one can see intuitively from Fig. [U(a) that 
all Wannier functions \W(k y ,x)) can be parameterized 
by one real parameter. If we define K y = k y + 2nx for 



k y e [0,27r), then 



\W Kv=hy+2 , x ) = \W{ky,x)) 



(5) 



is continuous in K y £ R. The center-of-mass position of 
jWOfj,) versus K y is shown in Fig. Q] (b). In this notion 
one can see clearly that the Wannier functions labeled by 
K y are analogous to the lowest Landau level wavefunc- 
tions of the ordinary QH state in the Landau gauge in the 
form of ^ Ky {x,y) = ( 7r /|i^-i/4 elX ^_ (x _^i|)V2/ 2 Bj 
which are also eigenstates of k y and localized in x di- 
rection. Actually the Wannier functions Wif B ) reduces 
exactly to the Landau level wavefunctions if we apply 
this formalism to the Hofstadter model|24| and take the 
limit of small magnetic field lg ^> 1. 

Although the definition of |Wif B ) seems to be simply 
relabeling the Wannier functions, it plays a key role in 
understanding the QH and QAH states on equal footing. 
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Once the analog of Landau level wavefunctions | Wk v ) is 
found, all many-body wavefunctions of FQH states can 
now find an analog in FQAH states. 

Laughlin states and pseudopotential Hamiltonians.- 
The FQAH analog of Laughlin states in FQH can be 
constructed by using the basis | Wr y ) to replace the Lan- 
dau gauge wavefunctions in the Landau level problem. 
The Laughlin wavefunction on a cylinder is given by [HI 



*i /m (2i) = nf 



,2-KZi/Ly 



D 27T2 i /L ! 



i<3 



(6) 



with fl a normalization factor, and Zj_ = x. L + iy^. If we 
define the wavefunction in occupation number basis as 



®({ n i}) = jw j Y[dxidyiip : 
v J i 



2-rmi/L v (•''J l/m( z i) (7) 



then the FQAH version of the Laughlin state can be de- 
fined as 



\i/m) = 52*({m})]l\w 3mi/Ly ) 



(8) 



{«<} 



It is straighforward to verify that such a state is invariant 
under the lattice translation symmetry. In the occupa- 
tion number basis, the wavefunction is the same as that of 
the Laughlin state, so that one can also define the FQAH 
version of the pseudo-potential Hamiltonian [3, Ell, for 
which the state \l/m) is a unique ground state. For 
example, for the v = 1/3 Laughlin state the pseudo- 
potential Hamiltonian can be written in the following 
second-quantized formjiil. [26j: 



H = UJ2 bib n 
nez 

with b n = 

ZeZ,n — I even 



(9) 



2 I t~ 



in which c n is the annihilation operator of the single 
particle state \Wk —2 wn /L y }- hi other words, c\ L \<j>) = 
\W2- K n/L y ) with \4>) the vacuum. The state |l/3) satisfies 
b n 1 1/3) = 0, which is thus the ground state of H (for 
U > 0). 

It is essential to show that the Hamiltonian © is in- 
deed a local Hamiltonian of the 2d lattice system. On this 
purpose it is convenient to consider the coherent states 



|z> = E e ~ 



■nz!-i 2i- 



\W. 



2-xn/L 



(10) 



with z = z\ + 1Z2 a complex variable and z\, Z2 the real 
and imaginary parts. The coherent state \z) is periodic 
in z — > z + iL y which shows that the variable z is defined 
on the cylinder. The coherent state is a superposition 
of K y eigenstates around K y = 2im/L y ~ 2-kzi, 



with the width of the distribution AK y ~ 2y / 7r. Since 
the x position is proportional to K y as x ~ K y /27r, one 
can see that the coherent state is local (in the sense of 
exponential decay) in both x and y directions. For the 
two-band Hamiltonian Q, the coherent state wavefunc- 
tion is shown in Fig. [T] (d) . 

The annihilation operator of the coherent state can be 
defined as 



c(z) 



nez 



- 7r ( 21 -^)' 



(11) 



Since the coherent state wavefunction is local, c{z) is a 
local operator which is a superposition of the real space 
annihilation operators Cj around the center-of-mass posi- 
tion (zi, Z2). The Hamiltonian (|9]) can be written in c(z) 
as 



H 



ULy 

4tt 2 



J dz 1 dz 2 b f (z) ■ b(z), b(z) = c(z)(-iV)c(z) 



Since c(z) and b(z) are local operators in the 2d lattice 
model, so is H. 

Construction of more generic wavefunctions. - The ap- 
proach discussed above can be easily generalized to ob- 
tain more general FQAH states, such as the Moore-Read 
state 27( with non-Abelian quasiparticles. In general, any 



FQH state that can be written in the occupation num- 
ber basis on the cylinder geometry has an analog in the 
FQAH system. If the FQH state is an exact ground state 
of some Hamiltonian, a similar Hamiltonian can be ob- 
tained for the FQAH state by using occupation number 
basis and/or coherent states. For example, in Ref. [28| 
a large class of FQH states have been constructed by 
using Jack polynomials. For such states the wavefunc- 
tion in occupation number basis is recursively known, 
so that the generalization to FQAH states can be done 
straightforwardly. Similarly, other approaches of obtain- 
ing wavefunctions in FQH system (e.g. Ref. [29j| ) can 
also be generalized to FQAH system. Our result pro- 
vides a systematical approach to obtain two-dimensional 
chiral topological states on a lattice model with full lat- 
tice translation symmetry and a small number of single 
particle states per unit cell. (On comparison, the ordi- 
nary FQH states can be considered as the FQAH states 
in a Hofstadtcr model in the continuum limit Ib 3> a 
with a the lattice constant and Ib the magnetic length. 
In this limit the number of states per magnetic unit cell 
is diverging cx l B /a 2 .) Our approach can also be fur- 
ther generalized to more generic QAH states with Chern 
number C\ > 1, and/or multiple occupied bands. 
Generalization to fractional topological insulators. - 
The approach of Wannier functions can also be further 
generalized beyond FQH states. In recent years, topo- 
logical insulators (TI) have been proposed and experi- 
mentally realized, which can be considered as generaliza- 
tion of QH states in time-reversal invariant systems. (30j 
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In two-dimensions, TI can also be generalized to frac- 
tional TI 3jJ, |32( , which are fractionalized topological 
states protected by time-reversal symmetry. The sim- 
plest model for fractional TI consists of two decoupled 
FQH states, formed by spin up and down electrons, with 
opposite Hall conductance. For example for Laughlin 
state the corresponding fractional TI state is \FTI) = 
1 1/m, f) (8) | — 1/m, \) . More generally the spin S z conser- 
vation can be broken by spin-orbit coupling, and the TI 
and fractional TI states remain stable as long as time- 
reversal symmetry is preserved. 32, 33| However, in the 
wavefunction approach it is difficult to construct a frac- 
tional TI wavefunction for a system without spin con- 
servation, since the wavefunction of spin up (down) elec- 
trons in the spin-conserved case are holomorphic (anti- 
holomorphic) functions, and it is difficult to introduce a 
translation invariant coupling between them. The Wan- 
nier function approach provides a natural way to con- 
struct a fractional TI wavefunction for a generic Hamil- 
tonian without spin conservation. 

On this purpose we generalize the expression of maxi- 
mally localized Wannicr functions in Eq. ^ to a system 
with N occupied bands (20j: 



(a) 10 



\W(ky,x)) 



1 



E 



p e -»/o** a, x (p m ,ky)d P:L 



\n,k x ,k y ) (12) 



with e %Bn ( k v\ i 
Wilson loop operator 



1,2, ..N the eigenvalues of the 



and 
the U(N) 



the corresponding cigenstates. a x is 
Berry's phase gauge field defined by a x (p x ,k y ) nm = 
—i(n,p x ,ky\-^-\m,p x ,k y }. Similar to the one-band 
case, the center-of-mass position of the Wannier function 
\W z (k y ,x)) is given by (x) = x - 6 l (k y )/2ir. For a TI 
with a nontrivial Z 2 invariant [33| in the occupied bands, 
the center-of-mass positions of the Wannier functions are 
doubly degenerate at time-reversal invariant wavevectors 
k y = 0, 7T, and there are two subgroups of Wannier func- 
tions with opposite odd winding number in the evolution 
of k y from to 27r.[l!| [2^ As an example, we consider 
the Bernevig- Hughes-Zhang (BHZ) model for HgTe topo- 
logical insulator I34f and also include the bulk- inversion- 
asymmetry term[35[ such that the model is sufficiently 
generic and has no spin conservation. The lattice Hamil- 
tonian can be written as 

h(k) = (M + 2B(2 - cosk x - cosfc y ))l ® r 2 

+ sin k x a z (g t x + sin k y l <g> r y + Adj, <S> T y (13) 

in which o~ a , r a are Pauli matrices in spin and orbital 
indices, and the time-reversal transformation is defined 
as T = ia y K with K the complex conjugation. The 
Wannier functions can be obtained for this model, as 
shown in Fig. [5] (a) . In the TI phase the two Wannicr 
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FIG. 2: (Colour online) (a) The center-of-mass position (x) 
of the Wannier functions versus k y , and (b) (x) versus the ex- 
tended wavevector K y (see text), for the topological insulator 
model dT3j with the parameters M = -1, B = 1/2, A = 0.2. 



functions wind around the k y circle to opposite direc- 
tions. Similar to Eq. ([5]), we can extend the defini- 
tion of k y by continuity, and define the Wannier func- 



1.2 



tions Wk") = I W 1,2 (k y , £1,2)), with K y = k y ± 2^X1,2 
for i — 1,2 respectively The center-of-mass of these 
relabeled Wannicr functions are shown in Fig. [S] (b). 
The important property of this Wannicr function basis 
is that the occupied states are decoupled to two bands 
of Wannier states, each of which has the winding prop- 
erty as the Wannier states of a QAH system (with oppo- 
site winding direction), even though the Hamiltonian does 
not preserve any spin conservation. By using the Wan- 



nier functions 



W 



1,2 

A" 



we can construct wavefunctions for 



fractional TI in the form of \FQAHi) ® \FQAH 2 ), which 
are direct product of FQAH state \FQAH\) formed by 

W\ \ and the time- reversed state \FQAH 2 ) formed by 

W\ ^ ■ Although it appears to be a direct product state 

of two FQAH states, in the physical spin basis, spin up 
and down electrons are entangled and the state cannot be 
written as a direct product state. When a time-reversal 
invariant perturbation is considered (such as changing 
the term A in the Hamiltonian (fT3"]) ). the effect on the 
ansatz wavefunction can be described since the single 
particle Wannier states W^ 2 \ depends on the pertur- 
bation. Thus we see that through the Wannier function 
approach one can obtain generic wavefunctions for frac- 
tional TPs without spin conservation. 
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